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\.  ABSTRACT 

Crack  propagation  tests  were  conducted  on  a composite  modified  double- 
base (CMDB)  propellant  with  the  use  of  center-cracked  strip  biaxial  specimens. 

Constant  strain  rate  tests  were  conducted  at  several  temperatures  (4.0°  to 
105°  F)  and  crosshead  rates  (0.02  to  200  in. /min)  to  define  the  crack  initia- 
tion and  propagation  characteristics  for  monotonically  increasing  strain 
history.  The  tests  were  conducted  at  ambient,  250,  and  500  psig  pressur'- 
to  evaluate  the  effect  of  pressure  on  initiation  and  crack  velocity.  A 
second  series  of  tests  was  conducted  to  evaluate  the  effect  of  a prestrain 
damage  history  on  crack  propagation.  In  the  second  series,  the  samples 
(without  precut  cracks)  were  initially  prestrained  to  15  to  25  percent  and 
held  for  a period  of  time  to  Induce  material  damage.  After  load  release 
and  sufficient  recover/  time,  cracks  were  inserted  in  the  specimens  and 
they  were  then  pulled  to  failure  at  a constant  strain  rate.  Similar  tests 
were  conducted  on  round,  notched  tensile  samples  to  define  the  critical 
stress  intensity  factor  (Kj^)  and  to  provide  a comparison  between  uniaxial 
and  biaxial  fraction  initiation.  Schapery’s  viscoelastic  fracture  theory 
was  used  to  evaluate  the  crack  velocity  data  under  constant  strain  rate 
conditions.  One  important  result  of  the  study  was  the  finding  that  tlie 
crack  velocity  depended  rather  strongly  on  imposed  strain  level. 

r: 

INTRODUCTION 

Crack  initiation  and  propagation  in  polymeric  materials  has  been  givvn 
considerable  attention  during  the  last  15  years,  stimulated  primarily  bv 
their  unique  viscoelastic  behavior  and  also  by  the  use  of  polymeries  as 
solid  propellant  rocket  fuels.  Filled  polymers  are  used  extensively  a;; 
solid  propellant  grains  which  must  undergo  a variety  of  environmental 
loading  conditions,  during  motor  storage  and  handling  and  during  actual  ^ 

operational  firing  conditions,  which  impose  pressure  loads  at  high  ratt.s 

and  for  relatively  long  times.  The  consequence  of  a crack  can  be  the  4 

failure  of  the  rocket  motor  as  a result  of  overpressurization,  case  burn-  i 

through,  erratic  pressure-time  response,  or  a variety  of  events  related  to 
structural  or  ballistic  performance  failure.  The  primary  work  over  tlie 
last  decade  has  been  aimed  at  defining  the  laws  which  govern  crack  initi- 
ation, propagation,  and  trajectory  under  motor-like  operational  ccud  1 1 ioij.'i 
to  arrive  at  an  assessment  of  overall  crack  criticality  for  a given  rocket 

motor  system.  | 
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A number  of  investigators  have  studied  viscoelastic  crack  propagation  1 
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along  the  lines  of  classical  elasticity  using  the  Griffith  solution  for 
crack  initiation.  Bennett  (11,12)  presented  the  first  known  data  on  solid 
propellant  crack  initiation  in  a state-of-the-art  composite  propellant  and 
noted  that  the  fracture  energy,  f,  was  time-  and  temperature-dependent. 

Use  of  the  data  for  motor  applications  was  made  on  the  basis  of  an  energy 
balance  using  finite  element  analysis.  Correlation  of  the  data  with  that 
from  analysis  using  the  thermodynamic  power  balance  and  the  Williams' 
spherical  flaw  model  was  not  quantitatively  satisfactory  (4,11).  No 
attempt  was  made  to  monitor  crack  velocity  during  the  tests.  Only  crack 
initiation  was  monitored  by  visual  observation. 

H Knauss  (5-8,13)  and  Schapery  (14,15)  have  made  significant  contribu- 

tions to  theoretical  portion  of  the  problem  of  viscoelastic  crack  propaga- 
tion. Knauss  (5,8)  developed  an  approximate  solution  for  a line  crack 
problem  which  incorporates  a failure  zone  length  to  introduce  time  and 
velocity  effects.  Knauss  and  Mueller  (5,6)  applied  the  theory  to  an 
I unfilled  polyurethane  polymer  with  encouraging  results  achieved  for  a 

, variety  of  loading  environments. 

i Using  a stress  intensity  approach  to  fracture  of  solid  propellant 

grains,  Francis  et  al  (16,17)  conducted  a program  on  an  unfilled  poly- 
urethane polymer  (Solithane  113),  an  epoxy,  and  later,  a PBAN  composite 
propellant  formulation.  The  primary  contributions  of  their  work  was  the 
discovery  that  superimposed  pressure  causes  about  a 30  to  40  percent 
increase  in  the  stress  intensity  factor  at  a given  crack  velocity,  thus 
I requiring  more  energy  to  propagate  the  crack.  This  fact  was  essentially 

known,  or  assumed  from  superimposed  pressure  tensile  tests,  but  was  not 
j directly  related  to  actual  crack  propagation  experimental  data  per  se. 

The  thermal  and  pressure  loading  data  on  small,  two-dimensional  samples 
of  Solithane  113,  obtained  by  a stress  intensity  factor  approach,  were 
encouraging. 

Swanson  (18)  was  the  first  to  apply  fracture  mechanics  to  composite 
modified  double-base  (CMDB)  propellant.  He  performed  a stress  analysis 
of  small,  precracked,  subscale  STV  motors  formulated  in  terms  of  stress 
intensity  factors.  Using  quasi-elastic  analysis  and  time-dependent 
critical  stress  intensity  factors,  he  was  reasonably  successful  in  pre- 
dicting crack  initiation  in  the  ST\^  subscale  motors  subjected  to  high-rate 
pressurization  tests. 

The  problem  of  crack  propagation  in  solid  propellants  has  onlv  heen 
briefly  mentioned  in  previous  literature,  leaving  vast  areas  untouched. 
This  is  particularly  true  for  the  solid  rocket  motor  grain  itself,  whic.'i 
must  undergo  a wide  range  of  environments  and  loading  conditions  not 
heretofore  studied  in  the  laboratory  either  experimentally  or  ,ina  1 vti  < a H v 
Schapery  (15,16)  has  rece  .'-ly  been  working  on  a mitre  genera  1 iced  appvoa.  i. 
which  may  be  able  to  resolve  uiore  complicated  load  histories  representa- 
tive of  actual  motor  conditions.  Schapery  generalized  the  Barenhlatr 
model  (1*^)  for  elastic  fracture  to  linear  visioelastic  conditions  and 
then  developed  the  governing  equations  for  crack  vi-Jocitv  in  terms  o!  t!,i 
linear  viscoelastic  material  properties. 
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Swanson  (20)  recently  used  the  Schapery  theory  to  evaluate  published 
data  on  a FBAN  composite  propellant  and  found  that  there  was  a very  good 
agreement  between  theory  and  experiment  over  a large  range  of  variables. 
Swanson  found  a time-dependent  fracture  energy,  T,  for  both  crack  initia- 
tion and  crack  propagation  must  be  assumed  whereas  P constant  could  be 
assumed  for  Sollthane  113. 

The  present  program  was  aimed  at  evaluating  the  Schapery  theory  when 
It  was  applied  to  a CMDB  propellant  under  similar  loading  conditions.  CMDB 
propellant  generally  Is  different  than  composite  propellant  as  it  has  a 
granular  filler  material  of  a characteristic  length  of  0.05  to  0.07  Inch 
and  uses  a nitrocellulose  binder  matrix  instead  of  polyurethane  and  poly- 
butadiene.  Aside  from  the  more  standard  loading  environments  of  tempera- 
ture, pressure,  and  rate  (time),  it  was  also  of  interest  to  compare 
multlaxlal  effects  (uniaxial  versus  biaxial)  and  investigate  the  effects 
of  prestraining  damage  on  crack  propagation  characteristics.  The  latter 
is  of  interest  to  predict  the  behavior  of  a crack  in  a motor  which  must 
undergo  previous  loading  history  before  actual  ignition  loads,  a situation 
which  is  the  case  for  most  all  missile  systems  in  use  today. 

THEORETICAL  BACKGROUND 


Schapery's  Theory  of  Viscoelastic  Fracture 

Schapery's  theory  of  viscoelastic  fracture,  as  well  as  some  applica- 
tion to  available  polymeric  materials.  Is  discussed  in  References  15 
and  16.  For  completeness,  the  main  highlights  are  reviewed  in  this  para- 
graph, since  the  CMDB  propellant  data  will  be  considered  with  respect  to 
the  existing  theory  in  later  paragraphs. 


Schapery  develops  the  theory  for  the  time-dependent  size  and  shape 
of  cracks  In  linearly  viscoelastic,  isotropic  media  for  Mode  I (tensile) 
opening  stress  conditions.  The  proposed  model  assumes  that  failure 
occurs  in  a small  region  behind  the  crack  tip  and  is  essentially  equiva- 
lent to  the  Barenblatt  model  (19),  except  that  no  restriction  is  made  on 
the  constitutive  properties  of  the  material  in  the  failure  zone.  Linear 
viscoelastic  theory  predicts  the  boundary  around  the  crack  tip  to  be 
cusp-shaped,  a geometry  which  is  based  on  the  assumption  of  linear  visco- 
elastic behavior  In  the  surrounding  material  and  which  is  demonstrated  in 
several  solid  propellants.  The  singularity  in  stress  at  the  crack  tip  due 
to  the  cohesive  forces  in  the  failure  rone  is  equated  to  the  negative 
of  the  singularity  in  stress  at  the  crack  tip  due  to  the  external  applied 
load,  so  that  the  resulting  stress  is  everywhere  finite.  Schapery  calcu- 
lates the  work  done  on  the  cohesive  (failure)  zone  by  the  surrounding 
linear  viscoelastic  material  and  equates  this  to  the  fracture  energy,  r. 
The  relationship  for  the  crack  velocity  is  then  developed  in  the  form 


I' 


(1) 


where : 


r = Fracture  energy  required  to  produce  one  unit  of  new  surface 
area 
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= Related  to  creep  compliance  of  surrounding  material 

t^  = Characteristic  time 

“ Opening  mode  stress  Intensity  factor 

The  characteristic  time,  t^,  is  related  to  the  length  of  the  failure 
zone , a , and  the  crack  velocity,  a,  by 


t 


a 


(2) 


where  the  time  may  be  viewed  as  that  required  for  the  crack  to  traverse 
the  failure  zone  at  the  crack  tip.  Cy  is  proportional  to  the  uniaxial 
plane  strain  creep  compliance  which,  for  a constant  value  of  Poisson's 
ratio,  is  given  by 


C^(t)  ^ 4(1  - W^)D(t)  (3) 

Compliance  (and  modulus)  data  for  polymeric  materials  may  often  be 
expressed  as  a power  law  in  time.  For  a one-  or  two- term  power  law  repre- 
sentation of  C,^,  equation  (1)  can  be  solved  explicitly  for  the  crack 
velocity  as 


where : 


P 2n  2n  23+n 
m 1 


1/n 


2(l+l/n) 


(4) 


C = C,  t 
V 1 


a^i  Ij  = Constants  or  material  parameters  defined  by  Schapery 
for  the  failure  zone 


If  one  can  assume  that  the  fracture  energy,  I',  is  a constant,  then  the 
propagation  law,  equation  (4),  mav  be  written  in  the  form 


a 


where  A becomes  the  term  in  the  brackets  In  equation  (4)  for  i:on,slaiii 

values  of  r and  I,  . * 

m 1 

Swanson  (20)  discusses  some  of  the  possibilities  oi  equatii.n  (5) 
if  r is  not  constant  and  evaluates  PBAN  propellant  data  for  a time- 
dependent  value  of  r.  However,  the  key  result  of  Schapery's  analysis  is 
that  crack  propagation  fracture  energy  cannot  depend  on  higher  order 
derivatives  of  crack  motion  or  histories  of  crack  motion.  For  purposes 
of  this  paper,  the  functional  dependence  of  I'  was  not  investigated  but 
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restricted  to  the  general  behavior  already  noted  in  the  original  objec- 
tives. Equation  (5)  says  that  the  relationship  between  the  stress  Intensity 
factor  (Ki)  and  crack  velocity  (a)  is  linear  when  plotted  on  log- log  paper 
and  has  a slope  proportional  to  the  exponential  (q) . 

Extension  of  Theory  to  Constant  Strain  Rate  Tests 

Since  the  data  to  be  presented  in  subsequent  sections  were  derived 
from  constant  strain  rate  tests,  the  crack  growth  law  will  be  reformulated 
for  the  particular  test  at  hand. 


In  the  constant  strain  rate  test,  the  applied  stress,  a,  is  derived 
from  the  Boltzmann  superposition  Integral  for  Isothermal  conditions  as 


<T(t) 


(6) 


which,  for  the  power  law  relaxation  modulus  given  by 


Vi<'>  * h' 


-n 


O) 


results  in 


1-n 


ff(t) 


(1  - n) 


(8) 


where ; 


= Strain  rate 

Ej^,  n = Power  law  material  constants 

It  should  be  noted  that  the  applied  stress  (<y)  is  the  far  field  stress 
which  would  be  present  for  an  infinite  strip.  The  necessary  correction 
factors  for  finite  strip  width  will  be  discussed  in  the  following  para- 
graphs . 

For  a centrally  cracked  plate  of  infinite  width,  the  stress  is  related 
to  the  stress  intensity  factor,  K^,  and  the  half-crack  length,  a , by 
the  expression 


The  crack  growth  may  be  expressed  in  terms  of  the  applied  stress  (<r) 
by  integrating  equation  (4),  using  the  chain  rule  to  change  the  independent 
variable  from  time  to  stress,  and  using  the  relationship  between  time  and 
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stress,  equation  (8),  to  get 


a 

^o 


i - C4(a^)'-/>^(aP 


where  equation  (9)  has  been  used  to  express  the  stress  intensity  factor 
in  terms  of  applied  stress,  and  where 


r - nMR^ 


1 

1-n 

1 - n\  1 

R.  E,  i P 


(11a) 


P = q + 


1-n 


(11b) 


q = experimental  value  from 
Equation  (5) 


(11c) 


(lid) 


= critical  stress 


(lie) 


Equation  (11)  holds  when  the  fracture  energy  r is  a constant,  the  assump- 
tion used  to  evaluate  the  CMDB  data. 

Finite  Geometry  Correction  Factors 


The  determination  of  stress  intensity  factors  from  experimental  data 
generated  by  finite  width  centrally-cracked  strip  biaxial  specimens  requires 
geometry  correction  factors.  Isida  (21)  generated  correction  factors  loi 
finite  width  and  height  for  the  strip  biaxial  sample  sub)ected  to  several 
types  of  loading  boundary  conditions.  The  factors  given  in  (21)  were 
verified  and  extended  over  a wider  range  of  variables  by  using  the 
TEXGAP  (22)  finite  element  program.  The  values  of  determined  over 
the  range  of  interest  showed  very  good  agreement  with  the  original  Isida 
values.  TEXGAP  and  Isida  values  were  usually  within  1 to  2 percent  oi 
each  other.  A subroutine  was  written  for  eventual  data  reduction  and 
consisted  of  the  correction  factors  for  finite  width  and  height,  allow- 
ing for  crack  growth  and  boundary  displacement. 

EXPERIMENTAL  PROGRAM 

Crack  initiation  and  propagation  tests  were  conducted  in  a series 
of  experiments  which  involved  both  uniaxial  and  biaxial  tinsile  lists. 

The  uniaxial  tests  were  used  to  generate  fracture  toughness  (Ki^. ) i'm 
uniaxial  stress  field  conditions.  The  biaxial  tests  were  designed  to 
provide  fracture  toughness  (Kj^^,)  values  as  a comp.arison,  but,  moreover, 
the  tests  were  primarily  used  to  furni.sh  track  velocity  data  os  a func- 
tion of  the  imposed  stress  intensity  factor  (Kj,)  for  a variety  of 
environmental  conditions. 


The  propellant  used  in  the  program  is  a (.MDB  formulation  des i Kuated 
EJC.  The  propellant  v?as  taken  from  an  8-year  old  Polaris  second  stage 
motor.  Typical  uniaxial  tensile  properties  are  195  psi  maximum  s'.ress, 

50  percent  strain  at  maximum  stress,  and  570  psi  initial  modulus.  The 
experimental  program  which  was  carried  out  is  discussed  in  the  following 
paragraphs . 

Uniaxial  FracU-re  Toughness 

The  sample  used  to  generate  uniaxial  fracture  toughness  values 

is  shown  in  Figure  1.  The  sample  is  rough-cut  oy  saw  from  hulk  propel- 
lant removed  from  a sectioned  motor  and  subsequently  machined  on  a lathe 
into  the  desired  configuration.  The  crack,  in  the  form  of  a circumferen- 
tial notch,  is  inserted  to  a known  depth  with  an  X-acto  knife  mounted  on 
the  lathe  fixture  during  the  machining  opcra*'iou.  Because  of  the  nature 
of  the  test  and  specimen  configuration,  there  is  n )t  hflicvt-d  to  ie  any 
significant  di  ' ference  in  as  a result  of  i !k  - .'•i  f ic  i,'  ' iricf  -,r. 

comparison  with  a natural  track.  Crack  growth  in  t lis  specimen  is  rela- 
tively last  and  should  not  indicate  any  difterem-;  it.  heh,i',;.>r  as  a 
result.  Once  the  samples  arc  machined  they  are  t .i* - ' ■ 'nocf  to  iixtai  tali 
cylinders  prior  co  testing. 


The  samples  were  pulled  at  crosshead  rates  ranging  from  0.C2  L,> 

200  inches  per  minute  with  an  Instron  test  machine  for  the  lower  rates 
and  a closed-loop  (tydraulic  machine  for  the  higher  rates.  The  effect  of 
temperature  was  determined  over  a range  from  20°  to  105  F.  Since  pressorc 
is  often  a strong  Influence  on  the  fracture  characteristics  of  composite 
rubber-base  propellants,  the  tests  were  also  conducted  at  superimpose  1 
pressures  of  0,  250,  oOO,  and  1000  psig  with  a pressv;re  chamber  around 
the  s-xnple.  Care  was  taken  to  ensure  Chat  the  sample  was  not  dair;Cged  jue 
to  handling  and  pretest  loads  during  grip  insertion. 

Biaxial  Fracture  Characteristics 


The  crack  initiation  and  propagation  behavior  in  a biaxial  sfrt'ss 
field  was  defined  through  Che  us<  of  centrally  cracked  strip  biaxi  .i 
specimens  as  shown  in  Figure  2.  as  with  the  uniaxial  samples,  the  str,. 
biaxial  specimens  were  prepared  from  bulk  propellant  by  a similar  pro- 
cedures. The  specimens  were  machined  to  the  final  configure  Lion , sltC' 
which  the  end  tabs  were  bonded  on.  TlTe  end  tabs  pro'/ide  the  sample 
loading  tran.sition  as  well  as  the  biaxial  stress  field  due  to  the  lat«-ral 
bonding  constraint.  The  end  tabs  were  mounted  in  a special  grip  which 
prevents  sample  rotation  during  the  test,  a problem  of  considerable 
significance  in  crack  propagation  testing. 

In  this  particular  study,  the  crack  length  (.2a^)  was  selected  as 
1.5  inches.  This  length  corresponds  to  the  width  of  commercial,  sin-.-.lc 
edge  razor  blades,  which  wero  used  to  induce  the  artificial  crackr.  in 
the  specimens.  The  cracks  were  m.adc*  Lri  the  specimens  before  testing 
after  careful  alignment  of  the  cutting  blade  over  the  center  of  the 
specimen,  equidistant  from  the  edges.  A transparent,  plastic  grid  is 
superimposed  over  the  crack  ares  and  taped  in  place  near  the  grips  to 
provide  a reference  base  during  the  crack  growth  period.  The  crack 


r 


pattern  was  made  more  distinct  I'v  placing  ^ sheet  nf  white  paper  on  the  i 

back  side  of  the  tan-colored  specimen  to  assist  in  the  final  photo  analysis. 

Crack  initiation  and  propagation  were  recorded  and  assessed  through  the  ^ 

use  of  slow  motion  (low  rate  tests)  and  high  speed  tmoderati'-to-high  rate 
tests)  photography  at  frame  rates  up  to  5000  frames  per  second. 

The  primary  tests  were  conducted  at  crosshead  rates  ranging  from  ‘ 

0.02  to  200  inches  per  minute  with  the  Instron  and  closed- loop  hydraulic 
test  machines.  The  tests  were  run  at  temperatures  of  40°,  77°,  and  105  ‘ F 
after  an  equivalent  dumny  sample  (same  thickness)  indicated  temperature 
equilibration.  Tlie  effect  of  superimposed  pressure  was  also  investi- 
gated at  77°  F for  pressure  levels  of  250  and  500  psig.  The  crack  growth  j 

characteristics  for  the  thermal  tests  were  observed  by  mounting  a plexi-  j 

glass  door  on  the  temperature  chamber,  while  those  for  the  pressure  tests  1 

were  observed  through  a 4-inch  thick  plexiglass  window  on  the  pressure  j 

chamber.  Optical  clarity  in  the  pressure  tests  was  particularly  good, 
since  the  sample  was  located  relatively  close  to  the  window  and  external 
lighting  cou]d  be  readily  used.  Care  was  taken  to  ensure  chat  saiiij>le 
surface  heating  was  avoided  by  using  the  lights  only  during  the  camera 
recording  bursts  and  by  keeping  the  burst  periods  short  (10  seconds  or 
less).  The  reference  time  was  determined  by  the  camera  system  (eleccroui. 
frequency  counter)  or  by  a stopwatch  located  in  the  viewing  area. 

A second  series  of  tests  was  conducted  to  evaluate  Che  effects  of  pre- 
strain  material  damage  on  crack  propagation  behavior.  In  this  series,  the 
strip  biaxial  samples  selected  for  the  study  were  pulled  to  predetermined 
strain  levels  of  15  and  25  percent  at  the  crosshead  rate  of  2 inches  per 
minute  and  then  held  for  30  minutes.  The  samples  were  then  unloaded  and 

allowed  to  relax  overnight  before  retesting.  The  cracks  were  inserted  _ 

after  the  recovery  period  and  just  before  the  retesting.  The  test.s  were 
conducted  at  ambient  pressure  at  temperature.s  of  40°,  77°,  and  105°  F and 
c ro,'.- she  .'id  rates  of  0.02,  2,  and  200  inches  per  minute  by  the  same  pro- 
cedures u.sed  tor  the  virgin  samples. 

Stress  rela>'.p.tion  modulus  data,  necessary  for  evaluation  of  crack 
growth  hel'.avior,  was  available  from  Reference  23.  Hiese  data  were  'ibtaiued 
on  round,  end-bonded  samples  similar  to  the  notched  fracture  ti'ughness 
sample  shown  in  Figure  1.  Tiae  strain  revel  used  for  these  dat.u  was 
tvpically  1/2  to  1 percent  strain.  The  data  are  shown  fitted  to  the  powo  ■ 

law  ‘‘or"i  ^ 

F,(t)  = (See  Eq  7) 

where  Ei  and  n are  material  constants  for  Che  particular  CMDB  propellant.  ^ 

The  power  law  I'it  works  well  over  the  seven  decades  of  Lime  shown  in 
Figure  3. 

DISCUSSION  OF  RESULTS 

Unia.'iial  and  Biaxial  Fracture  Tmigimess  ; 


The  critical  stress  intensity  tactor,  or  fracture  tougliness  (Kj^j,), 
determined  from  Che  uniaxial  round  notched  and  strip  biaxial  sample  is 


shown  In  Figure  4 as  a function  of  the  tlme-to- fal lute , tj..  Several 
conments  should  be  made  concerning  the  method  of  data  reduction.  The 
value  of  Kjj,  for  the  notched,  round  tensile  specimen  was  derived  from  a 
solution  by  Bueckner  (24)  and  given  by 

K,  = <r  VirD  F(d/D)  (12) 

ic  mn 

where ; 

D = 0>itside  diameter 

d = notched  section  diameter 

/T  = Net  section  stress  determined  at  maximum  load 
"mn 

F(d/D)  = Geometrical  correction  factor  shown  in  (24) 

The  photo  analysis  of  the  strip  biaxial  samples  was  used  to  provide  the 
value  of  irora  tfie  time  at  which  the  crack  t»>as  visually  observed  to 

crack.  The  computer  subroutine  automatically  input  the  correction  factors 
tor  finite  width  and  height.  In  summary,  the  uniaxial  was  determined 
from  the  maximum  load  whereas  the  biaxial  was  determined  from  visual 
observation  of  crack  growth. 

There  does  not  appear  to  be  any  difference  in  behavior  between 
uniaxia’  and  biaxial  crack  initiation  in  the  CMDB  propellant.  The  two 
methods  ot  determination  should  not  be  significantly  different,  since 
the  .lotched  round  tensile  sample  would  be  expected  to  exhibit  rather  fast 
crack  growth  from  initiation  to  sample  separation.  Pressure  also  does 
nvit  significantly  affect  the  time- to- failure  for  the  CMDB  propellant. 

The  functional  “^onn  of  Figure  4 fits 


BK 

Ic 


where ; 


B = 9.2  (lO"^^) 


Fffect  of  Strain  Rate 


The  evaluation  of  the  load-time  traces  and  the  film  records  by  using 
the  correction  factors  tor  finite  width  and  height  results  in  the  expcii- 
menta]  relationship  between  stress  intensity  factors  and  crack  velocitv 
as  shown  in  Figure  5.  The  solid  lines  represent  the  data  from  four  rep] ■ 
cate  testa  at  any  one  given  rate  and  typically  show  less  than  10  percent, 
variation.  Any  individual  curve,  for  a given  rate,  follows  the  form 
derived  by  Schapery  in  equation  (5)  wi  tli  the  values  given  by 


"I 


where  depend  on  the  particular  strain  rate  employed.  This  is  unusual 
since  all  previously  reported  data  on  composite  propellants  indicate  the 
value  of  the  exponent  to  be  on  the  order  of  6 to  7 and  a continuous,  smooth 
curve  rather  than  a dependence  on  strain  rate  as  shown.  However,  if  one 
cross  plots  several  values  of  Kj^  and  a for  a fixed  strain  level,  then  it 
becomes  readily  apparent  that  a continuous  function  Is  found,  but  that 
the  crack  propagation  exhibits  a strong  nonlinear  dependence  on  strain 
level.  Hie  crack  propagation  law  now  takes  on  the  form 


7.50 


(15] 


where  the  value  of  Ag  depends  on  strain  level.  Tne  exponent  agrees 
reasonably  well  with  experimental  data  on  composite  propellants  as  well 
as  the  value  given  in  equation  (13)  from  the  crack  initiation  data 

There  are  two  important  aspects  of  Figure  5.  Tlic  first  is  that  the 
data  obey  Schapery’s  power  law  crack  propagatl<?n  law,  equation  (5),  over 
several  decades  of  crack  velocity.  The  second  and  more  important  aspect 
is  that  then  is  a strong  effect  of  strain  level  on  the  crack  propagation. 
This  result  has  not  been  reported  in  the  literature  but  is  expected  to 
be  rather  significant  in  analyzing  actual  rocket  motor  situations. 

Schapery  (25)  has  indicated  that  this  result  is  not  totally  unexpected 
on  the  basis  of  nonlinear  viscoelastic  behavior  exhibited  in  composite 
prope llan ts . 

Closer  examination  of  equations  (4)  and  (5)  would  indicate  that  the 
strain  dependence  could  very  easily  enter  in  through  the  creep  cotr.pliance 
or  the  description  of  the  failure  zone.  Further  work  is  needed  in  this 
area  to  uefine  this  dependence. 

Temperature  and  Pressure  Effects 


The  effects  of  temperature  and  superimposed  pie.ssuro  on  crack  propnga-  j 

tion  are  shown  in  Figures  b through  8.  Hie  effect  of  temperature  follows  \ 

the  familiar  time -temperature  shifting  principles  shown  in  Figure  7.  The  1 

date  were  shifted  experimentally  by  using  a single  rate  and  temporarily  j 

ignoring  the  strain  effect  noted  earlier.  The  correct  procedure  would  j 

entail  a definition  of  the  crack  propagation  law  for  a given  strain  ievci , ] 

say  10  percent,  followed  by  the  sliifting  procedure.  Pressure  .shows  the  j 

same  effect  .as  lowering  the  temperature,  namely,  increasing  the  stress  1 

intensity  factor  at  a given  crack  velocity.  The  same  effects  are  seen  | 

in  composite  propellants  (16).  A pressure  enhancement  of  about  25  to  i 

30  percent  on  the  stress  intensity  factor  is  realized  with  the  addition  ^ 

of  500  psig. 

Prestrain  Damage  Effects 

Tlie  tests  on  strip  biaxial  samples  subjected  to  various  levels  of 
prcotrain  were  conducted  primarily  to  assess  the  difference  in  crack  i 

propagation  as  a result  of  prior  load  history,  namely,  a constant  strain  ] 

history  to  some  level  t^.  Figure  9 shows  the  result  of  the  tests.  At  ’ 

the  highest  strain  rate,  200  in. /min  crosshead  rate,  no  difference  is  ! 


i ! 


1 


seen  in  the  crack  propagation  behavior  after  the  prestraln  danage  cycle. 
Even  for  prestrain  levels  as  high  as  25  percent,  no  significant  differ- 
ence was  exhibited  in  the  response.  However,  as  one  goes  to  lower  strain 
rates,  possibly  0.02  in. /min  crosshead  rate,  the  effect  of  the  prestrain 
cycle  begins  to  affect  the  crack  propagation  response.  Tlie  effect  of  the 
prestrain  is  to  decrease  the  stress  intensity  factor  at  which  the  crack 
will  propagate  for  a fixed  velocity.  Referring  to  Figure  5,  where  the 
effect  of  strain  level  on  crack  propagation  was  first  noted,  it  can  be  seen 
that,  at  the  high  strain  rates,  the  propagation  takes  place  at  a higher 
strain  level  (e.g.,  14  to  22  percent).  Conversely,  at  the  lower  strain 
rates,  the  crack  propagation  begins  at  much  lower  strain  levels  (e.g., 

<8  percent).  I'herefore,  it  would  not  be  expected  to  see  much  difference 
in  behavior  at  the  high  strain  rates  for  prestrain  levels  on  the  order 
of  15  to  25  percent.  The  lower  strain  rates  should  exhibit  a different 
behavior  due  to  the  amount  of  microcracking  which  resulted  from  the  15  to 
25  percent  prestrain  cycle. 

Comparison  with  Schaperv's  Tlteorv 

The  cjn.st.-mt  strain  rate  data  shown  in  Figure  5 wa.s  evaluated  by 
using  hohaper'’'s  thiory  (15,16)  to  predict  the  crack  growth  from  measured 
material  properties  and  experimental  stress  values  after  correcting  for 
finite  eeometry  effects.  Using  equation  (10),  after  reformulating  to  the 
more  convenient  form 

Aa  = a - a (16) 

o 

the  data  were  used  t(.>  predict  the  expected  crack  growth,  .^a.  It  should 
be  remembered  that  there  is  some  latitude  in  empirical  fitting  of  the 
dat.'i  in  view  of  the  assumptions  which  can  be  made  regarding  f.  Howe\’er, 
in  viewing  the  data  and  theory  as  shown  in  Figures  10  through  14,  the 
qualitative  fii  oi  the  data  to  the  theory  appears  to  be  reasonably  good, 
alth.nigh  further  wt'rk  is  necessary  in  light  of  the  strain  effect. 

CONCLUSIONS 

The  crack  propagation  behavior  of  CMDB  propellants  appear  to  be 
qua li t a t i\ e ly  similar  to  that  previously  reported  for  unfilled  polymers 
(Solithane  113)  and  PBAN  composite  propellants.  Both  propellant  systems 
exhibit  strc>ng  temperature-  and  pressure-dependence  on  crack  growth. 
the  CMDB  propellant  studied,  Che  uniaxial  and  biaxial  fracture  toughness 
(Kj^j,)  are  in  good  agreement,  indicating  little,  if  any,  difference  in 
multlaxial  behavior  for  crack  initiation  for  opening  mode  tensile  condi- 
tions. Comparisons  with  Schapery's  theory  shows  that  the  CMDB  propellant 
obeys  the  power  law  crack  propagation  law,  equation  (5),  and  agrees 
qualitatively  for  the  stress-crack  growth  predictions  shown  in  Figures  1) 
through  14 . 

Perhaps  the  most  Important  results  discussed  in  this  paper  are  iht 
strong  effect  of  strain  on  crack  propagation  and  the  effect  of  prestrain 
damage  cycles  on  low  strain  rate  crack  propagation  behavior.  Both  of  these 
effects  are  extremely  important  if  one  is  to  adequately  characterize  the 


material  and  make  reliable  grain  predictions  under  thermal  or  pressure 
loading  conditions. 
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